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Abstract 

Spontaneous emission (SE) from a two-level atom in a photonic crystal (PC) with anisotropic one- 
band model is investigated using the fractional calculus. Analytically solving the kinetic equation 
in terms of the fractional exponential function, the dynamical discrepancy of SE between the 
anisotropic and isotropic systems is discussed on the basis of different photon density of states 
(DOS) and the existence of incoherent diffusion field that becomes even more clearly as the atomic 
transition frequency lies close to the band edge. With the same atom-field coupling strength and 
detuning in the forbidden gap, the photon-atom bound states in the isotropic system turn into the 
unbound ones in the anisotropic system that is consistent with the experimental observation in 
Phys. Rev. Lett. 96, 243902 (2006). Dynamics along different wavevectors with various curvatures 
of dispersion is also addressed with the changes of the photon DOS and the appearance of the 
diffusion fields. 
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I. INTRODUCTION 



Photonic crystals (PCs), a new class of optical materials with periodic dielectric structure, 
provide a way to control spontaneous emission (SE) through redistributing the photon den- 
sity of states (DOS) near photonic band gap (PBG). This control offers the key technology 
of manipulating light, such as in light emitting devices pj], quantum information processing 
[2I, or solar cells [3]. The special photon density of states near band edge changes the optical 
behavior of an atom in a photonic crystal including the appearance of photon-atom bound 
, spectral splitting [9|], enhanced quantum interference effects [lO|, co- 



states 
herent control of SE 



non-Markovian effects 



12|, etc. Among these studies, the 



photon reservoir of a photonic crystal could be well described by one band edge frequency 
ujc at a certain edge wavevector and a dispersion relation which leads to the special photon 
30S. The dispersion relation near the band edge is assumed to be isotropic in early studies 



12|. In a real three-dimensional photonic crystal with an allowed point-group 



symmetry, the photonic band structure is highly anisotropic, namely, the equal frequency 
surface near the band edge is no longer spherical. A vector form of photon dispersion re- 
lation is required to describe a more realistic picture of the band edge behavior. In this 



vector form of dispersion relation, the photon DOS is proportional to ^uj — uj^. while that 



of an isotropic dispersion relation is proportional to 1/ \/uJ — cUc, where u stands for the 
eigenmode frequency and Uc for the band edge frequency |13i]. This discrepancy of DOS will 
cause dynamical difference of spontaneous emission between the anisotropic and isotropic 

diffusion field 



systems. The most prominent difference is the existence o: 
was predicted theoretically by Yang et al 



16 



Ij,ll5|,ll6|,ll7|. It 



171] that the bare atomic transition 



frequency lying in the region near band edge will be shifted into the forbidden gap by the 
interaction with the radiation modes where an atom-photon bound state is generated. That 
is, the emitter with frequency lying in the forbidden gap will not yield SE. Recently, Barth 
et al. observed experimentally [l8|] that the anisotropic properties of a PC could be detected 
by employing single emitters such as quantum dots (QDs). When an emitter is placed in- 
side a PC with the anisotropic band structure, the additional anisotropy will imprint on the 
spontaneous emission of the system if the emission frequency lies in the forbidden gap. In 
this observation, the CdSe/ZnS quantum dots embedded inside artificial colloidal opals with 
direction-dependent band structure gave fluorescence image (SE) with an extra anisotropy 
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but did not emit light if embedded inside a PC with weak anisotropy of band structure for 
the emission frequency of the QDs lying in the forbidden gap. There exists inconsistency 
between the theoretical prediction and experimental observation because SE appears in the 
anisotropic PC system even for the emitter's frequency lying in the forbidden gap. These 
inconsistency arouses our attention for the physical properties of the anisotropic PC system. 
We found that the presence of SE from the QDs in the anisotropic PC results from the 
dynamical difference of SE from the isotropic (or weak anisotropic) PC. 

In this paper, we study the dynamics of SE from a two-level atom embedded in a PC 
with anisotropic one-band model (Fig. 1). Fractional calculus is applied to solve the non- 
Markovian dynamics of the anisotropic optical system with threshold-like DOS. We found 
that the dynamical discrepancy of SE between the anisotropic and isotropic systems is the 
consequence of the different DOS in these two systems and the existence of the diffusion field 
in the anisotropic system. With the same atom-field coupling strength and detuning in the 
forbidden gap, the photon-atom bound states near the band edge of the isotropic reservoir 
turn into the decaying states in anisotropic system. This alteration leads to the presence of 
SE in the anisotropic PC system which is consistent with the experimental observation by 
Barth et al. [18]. The new topic of how the curvature of the anisotropic dispersion relation 
affects the dynamical behavior of the system has also been investigated. The dynamical dif- 
ference between the systems with various curvatures of dispersion is observed and discussed 
on the basis of the curvature- dependent DOS and coupling strength. 

The paper is organized as follows. In Sec. II, the basic theory of the anisotropic system 
is given. In Sec. Ill, we solve the kinetic equation of the anisotropic system through 
fractional calculus and express the analytical solution in terms of the fractional exponential 
function. The dynamical behavior of the anisotropic system is compared with that of the 
isotropic system. And the influence of curvature of the anisotropic dispersion relation on 
the dynamical behavior is also discussed here. Finally, we summarize our results in Sec. IV. 

II. BASIC THEORY OF A TWO-LEVEL ATOM IN AN ANISOTROPIC ONE- 
BAND PHOTONIC CRYSTAL 

When the system of a two-level atom coupled to the field reservoir in a photonic crystal 
with anisotropic one-band model is considered, the Hamiltonian for this atom-field interact- 
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ing system can be written as 

H = hj2i(T22 + ^ huj^ata^ + ih'^g^{atai2 - (T2ia^), (1) 

k k 

where cTij = {i,j=l,2) are the atomic operators for a two-level atom with excited state 
1 2), ground state |1), and resonant transition frequency UJ21', and at are the annihilation 
and creation operators of the radiation field; uj^ is the radiation frequency of mode k in the 
reservoir, and the atom-field coupling constant gj: = ^^^21^ [__A^^j 2 . jg assumed to be 
independent of atomic position with the fixed atomic dipole moment 621 = d2iUd- Here V is 
the sample volume, is the polarization unit vector of the reservoir mode /c, and eo is the 
Coulomb constant. 

In the single photon sector, the wave function of the system has the form 

|^(t)) = i?(t)e--*|2,{0}) + X]Q(t)e-s*|l,{lJ) (2) 

k 

with initial condition 5(0) = 1 and C'g(O) = 0. Here B{t) labels the probability amplitude 
for the atom in its excited state |2) with an electromagnetic vacuum state and C^{t) for the 
atom in its ground state |1) with a single photon in mode k with frequency tug. We got 
the equations of motion for the amplitudes by projecting the time-dependent Schrodinger 
equation on the one-photon sector of the Hilbert space as 

^B{t) = -Y,g,CS)e-'^^' (3) 

k 

jC^kit) = 9uB{t)e^'''^' (4) 

with detuning frequency VL^ = ujj: — U2i. By substituting the time integration of equation 
(4) into equation (3), we have the time evolving equation of the excited-state probability 
amplitude 

jB{t) = - j'^ G{t - r)B{T)dr (5) 

with the memory kernel G{t - t) = J2k 4^''^'^'""^ ■ This memory kernel is related to 
the dispersion relation of photon reservoir. For the anisotropic photonic band gap (PEG) 
reservoir, the dispersion relation has a vector form and could be expressed by the effective- 
mass approximation as 1^ 

Uj^^ ujc + A(k - kc) (6) 



where A = fuoc/kl signifies different curvatures in different directions with scahng factor / 
whose value depends on the nature of the dispersion relation near the band edge ujc- The 
related memory kernel of the anisotropic system could be further expressed as 

G{t -t) = I p(a;)e-^(^---)(*-^)rfa; (7) 



with the curvature-dependent density of states 



m 



pM = \/^^e(^-^j. (8) 

where G(u; — Uc) is the Heaviside step function. This density of states has threshold- like 
behavior near band edge Uc which results in the different dynamical behavior of this atom- 
field interacting system from that of free space. By applying the complex Fresnel integral 
x'^~^e~°'^dx = T{p)/aP [2^, we could integrate this memory kernel as 

nif _-r\ - /3^^V /^^^ j[3^/4-A,(t-T)l /q^ 

with the detuning frequency Ac = UJ21 — of the atomic resonance frequency co'21 from the 
band edg e uJc-, coupling constant (S^^"^ = {uj2idl-^^ujl^'^) / {24TTeQhc^) , and t > r in the long time 



limit 2l|]. Substituting this memory kernel into the time evolving equation (5) and making 



the transformation B(t) = e^^^^Dit), we have the kinetic equation of this anisotropic system 
as 

J ^1/2 i7r/4 rt Tj/ \ 

_z,w + ,A.z,w^^L_y^_L)_,,. (10) 

III. DYNAMICS OF SPONTANEOUS EMISSION 



In this section, we will use the mathematical method of fractional calculus to solve the 
kinetic equation of the anisotropic system and discuss the dynamics of SE on the basis of 
the obtained analytical solution. 



When the right-hand-side term of the kinetic equation (10^ 
it as a Riemann-Liouville fractional differentiation operator 
That is, 



is considered, we could express 



22 



23 



2J] with order u = 1/2. 



dr = V{-ll2)^B(t) (11) 



/o (t-r)3/2 ^ ' 'dt^^ 

with Gamma function r(a;). The kinetic equation thus has a fractional differential form as 

d 2/?l/2pj7r/4 Jl/2 

-D{t) + ^AMt) + -^J^J^.m = 0. (12) 



In order to solve this fractional kinetic equation, we manipulated the fractional operators in- 
cluding the integral operator (d~^/dt~^) and the fractional differentiation operator d^^"^ / dt^/'^ . 
Mathematically, the adopted manipulation is not unique provided that one could justify the 
function arrived at is the solution of the original fractional differential equation. The first 
step of our manipulation yielded 



J-l 2/?V2pW4 j-i/2 

D(t) - D(0) + ^A.^Z^(t) + ^^3/. dt-y- D(t) = 0. 



(13) 



The second step gave 



^-1/2 



(14) 



Here the initial condition -D(O) =5(0) = 1 has been applied. The probability amplitude 
D(t) could be solved by performing the Laplace transform of the fractional derivative in- 



cluding 



L 



L 



3172^(0) 



_dt-y 



t=0 



and 



L[t- 



-1/21 



£(1/2) 



{lb) 
(16) 

(17) 



where b{s) is the Laplace transform of D{t). These procedures gave the Laplace transform 
of D{t) as 

^(^) = s + zA, + 2/5V2eW45i/2//3/2- (18) 
In order to express this equation as a sum of partial fractions, we need to find the roots of 
the indicial equation Y"^ + 2/?^/^e*'^/^y/ /^/^ + iAc = 0, where we have converted the variable 
s^/"^ into Y . There are two kinds of roots for this indicial equation: one with different roots 
Yi 7^ Y2 and the other with degenerate root Y^ = Y2. For the case of different roots, D{s) 
could be expressed as 



b{s) 



1 



1 



{^s-Y,^ 



1 



Y2)J (Fi-Fa; 



with 



73/2 



+ 



p 



(19) 



(20) 
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and 



73/2 \Jp 



(21) 



For the degenerate case, we have 



P 



Ac. The indicial equation becomes 



^1/2 

+ 2——-—Y + = 



p/2 



or 



^1/2 i7r/4\ 2 

The partial fractions of D{s) can thus be written as 

1 



D s) 



2 • 



(22) 



(23) 



(24) 



The dynamical solution of the probability amplitude D{t) could be obtained by applying 
the inverse Laplace transform on the partial-fractional forms of D{s) for the two cases of 
different roots and degenerate root. The applied formula of the inverse Laplace transform 
include 



{^/s - a)_ 



Et{ +aEt{Q,a^) 



(25) 



and 



L(yi-a)2j 



2atEt (^-i, + (1 + 2ah) Et (O, a^) + aEt Q, a^^ (26) 



CO (li)" 



-e"* being the two-parameter fractional exponential 



with Et{a,a) = r Er=o r(a+n+i) - dt- 
function of variable t, order a, and constant a [241]. The analytical solution for the fractional 
kinetic equation (Eq. (12)) of the anisotropic photonic crystal system is thus obtained. For 

1 



D{t) 



X 



1^2^,(1/2, Yl) - Y^E^{l/2, Yl) + Fie^i^* - F^e^^* . (27) 



For a If = Ac, 



-tEJ-. 

•2' 



J9/2 



^1/2 i7r/4 1 

m = -2^^^^tEti^,tP/f) - ^^^Eti-,t/3/f' 



p/2 



+ (1 + 2tt(3lf)e^^'/f' - r-^^t'l' 



TT. 



(28) 



Here we have applied the relation of the fractional exponential function for special values 
Et(-l/2, a) = aEt(l/2, a) + r^/^/v^ and Et(0, a) = e"*. 



This analytical solution, which has so far not been obtained, determines the dynamical 
behavior of the atomic excitation B{t) and the amplitude of the radiation field which could 
be obtained via B{t) in a standard way 16|, l25|]. Different kinds of the indicial roots Yi and 



I2 gave different dynamical behavior of the system which depend on the lying regions of the 
atomic transition frequency ( see Eq. (20) and (21)). As the atomic frequency lies in the 
region U21 < ujc + (3/ (Ac < P/ f^), the square of these roots is pure imaginary and the 



dynamical solution has some non-decaying terms. Near band edge {0021 = ujc, = Y?^ 



2 )i 



these non-decaying terms interfere each other severely which leads to the decaying solution. 
These non-decaying terms oscillate individually with time and form atom-photon bound 
states as the atomic frequency lies deeply in the forbidden gap {uj2i « oJc). On the 
other hand, when the atomic frequency is in the higher-energy region UJ21 > uJc + (3/ p 
(region of allowed band), the square of these roots is complex and the solution decays with 
time quickly. These two regions of atomic transition frequency given by the roots of the 



indicial equation are the discussing basis of the previous studies [Ij, [151, [161, Il7|. It was 
predicted in the previous studies that the bare atomic transition frequency lying in the 
region to'21 < ujc + will be shifted into the forbidden gap by the interaction with the 

radiation modes where a photon-atom bound state is generated. That is, there will not exist 
spontaneous emission in the anisotropic PC system if the atomic transition frequency lies 
in the region near band edge {uj2i = ujc)- This result is inconsistent with the experimental 
observation in Phys. Rev. Lett. 96, 243902 (2006) where SE appeared with an extra angular 
anisotropy in the anisotropic PC system as the emission frequency of the embedded QDs 
lies in the forbidden gap. 

Instead of analyzing the integration contours for the probability amplitudes in previous 
solving procedures, we plot the dynamical behavior of this anisotropic system directly from 
the analytical solution. Based on the excited-state probability amplitude P{t) = \B{t)f = 
\D{t)f, the dynamical behavior of the anisotropic system is shown in Fig. 2. It could be 
observed from Fig. 2 that typical characteristic of non-Markovian dynamics including non- 
exponential decay and atom-photon bound states exists in the system which results from the 
special (threshold-like) density of states (see Eq. (8)). When the atomic transition frequency 
lies in the bandgap (A^ < 0), the system exhibits photon-atom bound states and decaying 
states in the allowed band (Ac > 0). The dynamical behavior of the anisotropic system 



is almost the same as that of the isotropic system [8| except for the smaller probability of 
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bound states and faster decaying behavior of decaying states. The dynamical difference of 
spontaneous emission in the two systems results from the different DOS in the two systems 
and the existence of diffusion field in the anisotropic system. As the atomic transition 
frequency moves from the bandgap to the allowed band, the density of states "seen" by 
the emitted photon is singularly large near band edge in the isotropic system and small in 
the anisotropic case. The singularity of density of states in isotropic system leads to the 
appearance of coherent propagating field while not large enough density of states results 
in the coexistence of incoherent diffusion field and coherent propagating field in anisotropic 
system. The energy transfer from the locahzed field to the diffusion field for the bound states 
of the anisotropic system leads to the smaller probability in excited level and coexisting 
energy of diffusion field and propagating field for decaying states results in faster decaying 
of the excited population. 

This dynamical difference in the isotropic and anisotropic systems is more obvious as the 
atomic transition frequency lying close to the edge of the PBG reservoir which is shown 
in Fig. 3. As the atomic transition frequency lying close to the band edge, states in the 
isotropic system exhibit bound (Ac < 0) or slow decaying (Ac > 0) behavior. In the 
anisotropic system, however, almost all of these states display fast decaying behavior. That 
is, the bound states close to the band edge of the isotropic system change to the decaying 
states which lead to the appearance of SE in the anisotropic system. This result is consistent 
with the experimental observation in Phys. Rev. Lett. 96, 243902 (2006). In order to 
investigate the local optical properties of PCs, Barth et al. doped the PCs of artificial 
colloidal opals with CdSe/ZnS core-shell quantum dots whose emission frequency hes inside 
the forbidden gap and the linewidth was narrower than the width of the band gap. They 
demonstrated that the characteristic patterns of fiuorescence image from different quantum 
dots carried information on the modification of the optical mode density which arose from 
the direction-dependent photonic stop band. The anisotropic band structure of the artificial 
colloidal opals, which corresponds to the direction-dependent photonic stop band, brought 
an extra angular anisotropy of fiuorescence image that was detected by defocused wide-field 
image of the single CdSe/ZnS quantum dots. These quantum dots did not emit light if they 
were embedded inside a PC with weak anisotropy of band structure. That is, spontaneous 
emission appears only in the anisotropic photonic crystal as the emission frequency of the 
embedded quantum dots lies in the forbidden gap, but they did not emit light with bound 



9 



states in the weak anisotropic PC system. The correctness of our results is vahdated by 
the presence of spontaneous emission in the anisotropic photonic crystal which differs from 
the prediction of the previous studies Q, Q, Q. Dynamical difference of the isotropic 
and anisotropic PC systems leads to the appearance of fluorescence image in the anisotropic 
system. Close to the band edge, DOS in the anisotropic case is nearly zero so that the 
incoherent diffusion field emerges to release some of the radiation energy. This diffusion field 
increases as the atomic frequency shifts from the bandgap to the band edge and decreases 
as the frequency shifts to the allowed band. The transferred radiation energy from localized 
field and propagating field reaches maximum at the band edge. The phenomenon of the 
dynamical behavior of the two systems arriving at the great difference close to the band 
edge illustrates the existence of the diffusion field in the anisotropic system. 

In Fig. 4, we show how the curvature of the dispersion relation affects the dynamical 
behavior of the anisotropic system which has so far not been explored. The solid lines 
are plotted for the system with the larger curvature (/ = 1) of dispersion relation and 
dotted lines for the system with the smaller curvature (/ = 0.8). For the bound states 
(Ac < 0), the excited-state probability P(t) of the system with the smaller curvature has 
the smaller value than that of the system with the larger curvature. On the other hand, the 
excited-state probability of the system with smaller curvature decays faster than that of the 
larger-curvature system for the decaying states (Ac > 0). Without changing the units of 
energy (coupling constant (3) and time (1//5), the dispersion relation of the smaller curvature 
has the larger DOS and coupling strength. As the atomic transition frequency moves from 
the bandgap to the allowed band, this larger DOS leads to the earlier appearance of the 
diffusion field that corresponds to the earlier energy transfer from the localized field to the 
diffusion field for the bound states and earlier coexisting energy of the diffusion field and 
the propagating field for the decaying states. This earlier energy transfer and coexistence 
resulting in the smaller value of probability in the bound states and fast decaying in the 
decaying states cause different optical behavior of the anisotropic system along different 
wavevectors. 
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IV. CONCLUSION 



We have used fractional calculus to solve the non-Markovian dynamics of the optical 
system consisting of a two-level atom coupled to a PBG reservoir with anisotropic one-band 
model. It is the first time in the anisotropic system that the analytical solution of the 
kinetic equation is obtained in terms of fractional exponential functions. The dynamical 
spontaneous emission has almost the same behavior as that of the isotropic system except 
for the smaller probability of the photon-atom bound states (Ac < 0) and the faster decaying 
rates for the decaying states (Ac > 0). The dynamical difference originates from the different 
DOS in the two systems and the existence of diffusion field in the anisotropic system. Not 
large enough DOS near band edge in the anisotropic system leads to the appearance of the 
incoherent diffusion field and energy transfer from the localized field to the diffusion field. 
The dynamical difference between the anisotropic and isotropic systems manifests itself more 
clearly as the atomic transition frequency lies close to the edge of the PBG reservoir. With 
the same atom-field coupling strength and detuning frequency in the forbidden gap, the 
bound states in the isotropic system turn into the unbound states in the anisotropic system. 
This change leads to the presence of SE in the anisotropic system that agrees with the 
experimental observation in Ref. Il8| where spontaneous emission happens only in the strong 
anisotropic PCs but not in the weak anisotropic system for the emission frequency lying 
in the forbidden gap. The presence of spontaneous emission in the anisotropic photonic 
crystal validates the correctness of our results while illustrates the inconsistency with the 

n n n 

prediction of the previous studies [IJ, Il5|, [IGj. The existence of the diffusion field in the 



anisotropic system is elucidated by the dynamical behavior of the two systems arriving at 
the great difference close to the band edge. We also investigated the new topic of how 
the curvature of the anisotropic dispersion relation affects the dynamical behavior of the 
anisotropic system. Without changing the units of energy and time, the dispersion relation 
of the smaller curvature has larger DOS and coupling strength which leads to the earlier 
appearance of the diffusion field and energy transfer. This earlier energy transfer resulting 
in the smaller and faster-decaying probability causes the different optical behavior along 
different wavevectors in the anisotropic PC system. 
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FIG. 1: (a) A two-level atom with excited state |2) and ground state The transition frequency 
i02i is nearly resonant with the frequency range of the PBG reservoir, (b) Directional dependent 
dispersion relation near band edge with edge frequency ujc ■ (c) DOS of the anisotropic one-band 
effective mass model. 
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FIG. 2: Dynamics of SE of the anisotropic PC system as a function of (3t for various values of the 
atomic detuning frequency (Ac = — Wc). States of Ac//? < inside the band gap and states of 
Ac/ /3 > within the allowed band. 
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FIG. 3: Dynamics of SE of the anisotropic (Aniso., solid lines) and isotropic (Iso., dashed lines) 
systems with detuning frequency close to the band edge of PBG reservoir (Ac = UJ21 — ^ 0). 
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FIG. 4: Dynamics of SE of the anisotropic system with two kinds of curvatures in the dispersion 
relation. Solid lines for the system with smaller-curvature dispersion relation and dashed lines for 
the larger-curvature system. 



16 



1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 





VP = -5 




A,/P = -1 




A^/p = 0.5 




a}p = 2 






— 


. 1^ 1 1 1 y T 1 '\->- r !■ T 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 

pt 



1.0 
0.9 
0.8 
0.7 

g-o.s 

0.4 
0.3 
0.2 
0.1 
0.0 



1 
1 

II 


A^p = - 0.1, ISO. 


II 
II 
II 
1 1 


A^p - 0.5, ISO. 


1 1 
1 1 
1 1 


A^p = - 0.1, Anise. 


1 1 
1 1 
1 1 

I 1 -.-H 

II / 


A^p - 0.5, Aniso. 






1 1 1 
1 1 / 
1 1 1 
1 \ / 
1 \ * 

1 

1 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 



1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 



: 


VP = 


-5, larger curvature 




VP= 


-0.5, larger curvature 




VP= 


2, larger curvature 




VP = 


-5, smaller curvature 






-0.5, smaller curvature 




VP = 


2, smaller curvature 




u 

- \\ 






- \ \ 






1 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



1 2 3 4 5 6 7 8 9 10 1 1 12 13 14 15 

Pt 



